Statistical thermodynamics of a two dimensional relativistic gas 
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In this article we study a fully relativistic model of a two dimensional hard-disk gas. This 
model avoids the general problems associated with relativistic particle collisions and is therefore 
an ideal system to study relativistic effects in statistical thermodynamics. We study this model us- 
ing molecular-dynamics simulation, concentrating on the velocity distribution functions. We obtain 
results for x and y components of velocity in the rest frame (F) as well as the moving frame (F'). 
Our results confirm that Jiittner distribution is the correct generalization of Maxwell-Boltzmann 
distribution. We obtain the same "temperature" parameter P for both frames consistent with a 
recent study of a limited one-dimensional model. We also address the controversial topic of tem- 
perature transformation. We show that while local thermal equilibrium holds in the moving frame, 
relying on statistical methods such as distribution functions or equipartition theorem are ultimately 
inconclusive in deciding on a correct temperature transformation law (if any). 
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I. INTRODUCTION 



The question of how thermodynamic properties trans- 
form in a moving coordinate system were raised soon af- 
ter Einstein's fundamental paper in 1905 In no more 
than half a century the introduction of several relativis- 
tically consistent generalization of thermodynamics led 
to such a confusing atmosphere in which one could not 
decide whether a moving body appears cooler, hotter, or 
at the same temperature as the body at rest. The most 
cited view is presented by Planck Q and Einstein 3, 
who beUeved that temperature of a moving body would 
be Lorentz contracted. A different view was proposed 
later by some authors notably Ott Q and Arzelies 
suggesting that a body in motion would appear rela- 
tively hot. Finally, in 1966 Landsberg 0, Q put forth 
the third suggestion, namely, the Lorentz-invariant tem- 
perature view. However, 30 years later Landsberg and 
Matsas d, [§| and recently Sewell [l3] proposed another 
view, that of nonexistence of universal Lorentz transfor- 
mation of temperature that further intensified the con- 
troversies over the subject. 

Since its early days relativistic thermodynamics has 
changed from a theoretically interesting problem to a 
practically important subject due to its application in 
the proper interpretation of experiments in high energy 
and astrophysics [O, [13, • Nevertheless, there is still 
no consensus on many features of this theory. One rea- 
son for the ongoing discussion is the lack of experimental 
evidences or numerical investigations. Among few ex- 
ceptions is an interesting paper by Cubero et al. 
who have shown that a simple one-dimensional model of 
relativistic dynamics favors Jiittner distribution function 
as the correct generalization of Maxwell-Boltzmann 
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(MB) distribution, 

/,7(v) = mS(v)2+'* exp[-/3jm7(v)]/Z;, (1) 

where d is dimension, Zj is normalization constant, E = 
m"f{-v) is relativistic energy, and 7(v) = (1 — w^)^^/^ is 
the Lorentz factor in natural units with speed of light 
c = 1. 

Although the model used in [3| is one-dimensional 
and lacks many features of a real physical system, it pro- 
vides strong evidence against other generalizations of the 
Maxwellian, especially the "modified" Jiittner function 

/m,/(v) - „ ; , exph/3M./m7(v)]. (2) 

Zmj m7(v) 

Jiittner distribution can be used as the cornerstone of 
our understanding of relativistic statistical mechanics in 
the same manner that MB distribution illuminates the 
underlying microscopic roots of classical thermodynam- 
ics. The most challenging step, however, is defining a 
proper thermometer in order to relate the Lagrange mul- 
tiplier P to the temperature of the system. This prob- 
lem is mostly treated as trivial in the literature but one 
should note that, the correct transformation of temper- 
ature, like any other quantity, depends crucially on the 
practical methods we implement for its measurement. 

Here, we model a two-dimensional (2D) gaseous sys- 
tem with realistic features which at the same time allows 
for implementation of full relativistic dynamics. Since 
this model is both realistic and fully relativistic, it can 
be used as an ideal numerical laboratory in order to 
investigate many issues concerning relativistic general- 
ization of statistical thermodynamics. Using standard 
relativistic transformations, we obtain directional distri- 
bution functions for both the rest as well as the mov- 
ing frame. We study these functions numerically using 
molecular-dynamics simulations of our 2D model. Our 
results indicate that Jiittner distribution is the correct 
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generalization. We also show that the same temperature 
parameter is obtained in both frames. Finally we discuss 
the implication of our results for a proper temperature 
transformation. In this regard, while verifying the impor- 
tant concept of local thermal equilibrium, we argue that 
these methods are ultimately inconclusive on deciding a 
correct temperature transformation law. 



II. MODEL 

We propose to study an idealized two-dimensional sys- 
tem of impenetrable hard disks with purely repulsive bi- 
nary interaction U (r) , 



U{r) = 



'CO, r < a 



0, 



r > a. 



(3) 



The disk-like particles move in straight lines at constant 
speed and change their momenta instantaneously when 
they touch at distance a [l^. Hence, in order to simu- 
late the dynamics, we must find the next collision and 
compute the changes in momenta of the colliding pair, 
considering the relativistic laws of conservation of energy 
and momentum in two-dimensional space. In order to 
solve this problem exactly, we add the assumption that 
when two hard disks collide, the force is exerted along 
the line connecting their centers, = — Yj. There- 
fore, the components of momenta perpendicular to 
remain unchanged — Pi,±) and the parallel com- 

ponents change in the same way as the one-dimensional 
case |20l, 
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FIG. 1: Equilibrium velocity distributions in the rest frame F: 
numerically obtained x component of single-particle velocity 
distribution(-l-) from simulation of A'' = 100 particles of mass 
m — 0.1. (a) Here, e = 2.28m and the corresponding temper- 
ature parameters are /3.j = 11.4, I3mj = 7.8. (b) e = 1.04m, 
l3j = 253.1 and Pmb = 259.8. A significant deviation from 
modified Juttner function is evident in the relativistic regime. 
Similar results are obtained for f{vy). 



III. RESULTS 



Rest frame 



E, = 7(«c™)'[(l + vl^)E, - 2w«,||], (4) 

where hatted quantities refer to momenta after collision 
and Vcm = {p i,\\~^Pj.\\) I {Ei~\-Ej) is the collision invariant, 
relativistic center-of-mass velocity of the two particles. 
With the same rules for particle j, a deterministic, time- 
reversible canonical transformation at each collision is 
defined. The additional assumption means that particles 
do not slide on each other when they collide. In contrast 
to the one-dimensional model, such elastic binary colli- 
sions lead to equilibrium even if colliding particles carry 
the same rest masses [14]. In our simulation we have used 
N particles of equal rest masses m that are constrained 
to move in a square box of linear size L. We use peri- 
odic boundary condition. Note that in order to simulate 
a stationary system in the rest frame, the center-of-mass 
momentum must be put to zero manually. This condition 
would automatically be satisfied (if not at each instant 
but at least on time average) if fixed reflecting walls were 
used [211]. 



In order to obtain the equilibrium state of the system 
we let the two-dimensional gas equilibrate (typically af- 
ter 10^ collisions) and measure velocities of particles at 
equal times with respect to laboratory frame. To collect 
more data, we repeated this procedure every \QN colli- 
sions. Simulation results for N = 100 particles are pre- 
sented. Particles are initially placed on a square lattice 
of constant L / \/N and velocities are chosen randomly in 
(5- vicinity (with 5 a small number) of |vo| — \/\ — (1/77)^, 
corresponding to mean energy per particle e — rim. As 
77 — !■ 1, the kinetic energy becomes small compared to the 
rest mass energy, recovering the dynamics of (classical) 
non-relativistic model. The 77 — s- 00 limit, on the other 
hand, corresponds to a highly relativistic model. 

Theoretically, the x(y) component of Jiittner and mod- 
ified Jiittner velocity distribution is obtained by integrat- 
ing Eqs.lH) and over Vyivr^^, 

2 

/j(«.) = ^7(«x)'[i^2(/3jm7(«x)) +i^o(/9,/m7(z;,))], 

(5) 
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fMj{Vx) 



2m 



Z 



MJ 



^{vxf Ki{(3Mjm^{vx)), 



(6) 



with Kn denoting modified Bessel functions of the sec- 
ond kind [i^. Here, the parameter I3j/mj is determined 
by means of the following procedure: we have, for the 
average energy. 



Etot/N : 



(i''v/(v)TO7(v). 



(7) 



|v|<l 



Computing the right hand side (rhs) of Eq.Q for Jiittner 
and modified Jiittner in the two-dimensional case gives 
rhsj = + 2(3m + 2)/(/3(/3m -t- 1)) and rhsMj = 

1//3 4- m, respectively. By inserting Etoti N and m 
into these equations, the parameter I3j/mj consistent 
with Jiittner and modified Jiittner velocity distribution 
is uniquely determined. 

We now check these results by considering two cases 
with r] — 1.04 and rj = 2.28. As shown in Fig.l, 
the obtained single particle distribution of velocity x- 
component agrees with Jiittner function (solid lines) 
in both regimes and converges to MB distribution in the 
non-relativistic limit. A significant deviation from mod- 
ified Jiittner distribution (dashed lines) is also evident. 
Exact same diagrams are obtained for y-component of 
velocity (not shown). In particular, the y-component re- 
sults were fitted with the same parameter /3 as the x- 
component data. This shows that our system has equili- 
brated properly through successive collisions. 



B. Moving frame 





FIG. 2: Equilibrium velocity distributions in the moving 
frame F' with relative velocity u — 0.5. The system pa- 
rameters are the same as Fig. 1(a) in particular f3j = 11.4 
and Pmj = 7.8. Note the breaking of symmetry in part (a), 
where the inset shows more details of the peak. 



and dashed lines are the velocity x{y) component of ve- 
locity distribution obtained, respectively, by integrating 
Eqs.lH]) and jSJ over Vy{vx),e.g., 



We now turn to the more interesting question of equi- 
librium velocity distribution of a relativistic gas in mo- 
tion. For this, we examine the system from the point of 
view of an observer who sees that the system's frame, F', 
is moving with a uniform velocity u in x-direction with 
respect to his rest frame, F. Using the entropy max- 
imization principle, the single-particle distribution will 
be determined by an additional constraint on the sys- 
tem, namely, that of a definite total momentum p' [23| , 

d f /\d+2 

f'A^') = exp[/3j7(zi)m7(v')(l - u.v')] (8) 



m''7(v')'^+^ exp[/3MJ7(")"i7(vO(l - "-vQ] 
^{u)Zmj 7(«)ni7(v')(l-u.v') 

(9) 

The primed quantities are measured in the moving frame 
and the additional 7(u) term in denominator is due to the 
contraction of the moving box that encloses the system 
[131. Figure 2 shows the results for a system similar to 
Fig. 1(a) with u = 0.5. Note that here, x{y) component 
of velocities are measured F'-simultaneously. The solid 



f'AO = — 7K)'[i^2(/3m7(u)7«)(l - <)) + 
Zj 

ifo(/3m7(«)7«)(l-0)],(10) 

Other components of these distributions cannot be ob- 
tained in closed form and are therefore plotted in Fig. 2 
using numerical integration. The parameter (3 used to fit 
data in Fig. 2 is the same as that of the rest frame. Note 
that this parameter is obtained exactly as a function of 
system parameters (to, e) in F. However, in F', /3 is a 
fitting parameter which turns out to be the same as that 
in F. 

Therefore as clearly seen from numerical results, our 
two dimensional model shows Jiittner distribution as the 
correct relativistic version of MB distribution. The fact 
that the same parameter /3 is obtained from both x and y 
component velocities shows equilibration. However, more 
importantly, the fact that same (3 is obtained from both F 
and F' frames seems to indicate the invariance of temper- 
ature consistent with earlier work of Landsberg Q and 
previous simulation results [l^l- We now discuss if this 
agreement can shed light on the long-lasting question of 
how temperature transforms in a moving frame. 
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FIG. 3; Local thermal equilibrium in the box as seen by 
a moving observer: The profile C{X,Y) is measured F'- 
simultaneously for the same system as in Fig. 2 with Lx and 
Ly, each divided into ten parts. The expected value of 
[/37(m)]-^ is 0.076 which agrees well with the obtained value 
of C throughout the lattice. 



there exists no universal and continuous Lorentz trans- 
formation of temperature at all [1, [§] . The argument is 
based on the fact that black body radiation of a moving 
body is direction dependent. Therefore, a bath which 
is thermal in an incrtial frame is non-thermal in a mov- 
ing frame. Does local thermal equilibrium (LTE) hold in 
our system? To check, we divide the box into n cells 
of area AA' = {L'^L'y)/ii? and calculate the quantity 

C =<C v'ixl'^'i ~ '^P'ix which we consider to be propor- 
tional to temperature, in each cell. The numerical result 
shown in Fig. 3, indicates that each cell, as seen by a mov- 
ing observer, is characterized by a constant value which 
coincides with the value of [/37(m)]~^. This indicates that 
LTE, which is the necessary condition to introduce a well 
defined temperature, is fulfilled at least for our model. 



IV. PROPER THERMOMETER 



CONCLUDING REMARKS 



A commonly used definition of equilibrium tempera- 
ture in the literature is T — (A;b/3)~^, where (5 is the 
Lagrange multiplier emerging in the velocity distribution 
function. One may use this definition and the equality of 
parameter /3 in moving and rest frame (Fig. 2) to deduce 
that temperature is Lorentz invariant, T' — T 14]. How- 
ever, there is another point of view which is also consis- 
tent with our results. Comparing Eq.® with the general 
form of the distribution function as / oc ^^i'^+P'^+'^-'P) ^ 
one is led to believe that the Lagrange multiplier is ^{u)j3 
which gives T' = T/j{u) 0,1, [2|. 

To see this point better, note that the above def- 
inition of temperature has its roots in the equiparti- 
tion theorem as well as basic thermodynamic relation, 
dE = TdS — PdV + lidN . The relativistic version of these 
methods is widely used in order to find the temperature 
transformation law 

0, i, S, Hi, m, ii] . However, despite 
the popularity of these approaches, they are not decisive 
either |26l. l27l. [28l| . To illuminate, consider the Lorentz- 
invariant equipartition theorem for a system moving with 
velocity u parallel to x axis , 



»'2 1 

Hit I "ill ^ 



(11) 



where the primed quantities are measured in the moving 
frame [e.g., m[ = 7(vi)mi] and averages ^ ... ^ are 
taken F'-simultaneously. One may apply either hypothe- 
sis that < pg/m'^ - up[^ > or < -i{u){pf^lm[ - up[^) > 
are the statistical thermometer of the moving system and 
find it compatible with the generalized theorem. Thus, 
as Landsberg has mentioned in "the argument from 
equipartition does not enable one to discriminate on the- 
oretical grounds between accepted theory and Lorentz- 
invariant temperature." 

Furthermore, in some recent papers it is claimed that 



It seems that the longstanding issue of relativistic ther- 
modynamics is related to the longstanding issue of irre- 
versibility in thermodynamics. The lack of consensus on 
these issues is related to the lack of concise mapping be- 
tween dynamical description of a system on one hand and 
a thermodynamic description on the other. We cannot 
define temperature (or entropy) as an exact function of 
dynamical variables. In this work we have modeled a 
useful, realistic system of a relativistic gas which over- 
comes the difficulties associated with implementation of 
particle interactions in a relativistically consistent man- 
ner [3, [2^ . We have shown that Jiittner function is the 
correct velocity distribution function in both rest and 
moving frames, with components either along or perpen- 
dicular to the relative velocity u. Furthermore, our re- 
sults indicate that, with a certain definition of statistical 
thermometer, one can choose (3' = /3, i.e., a Lorentz- 
invariant temperature, without running into inconsisten- 
cies. However, fi' = j{u)f3 could just as well be argued 
to be a valid choice, depending on a choice of thermome- 
ter. Such inconclusiveness inherit in statistical analysis 
like ours leads one to consider a covariant formulation of 
thermodynamics where temperature is generalized to a 
tensorial quantity whose transformation is no longer an 
issue [26l . [SO] . In this view thermodynamic temperature 
is considered as a proper feature of a thermodynamic sys- 
tem, much like mass in relativistic mechanics [SI]- 
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